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ABSTRACT 

We present the analysis of serendipitous sources in a deep, 500 ksec, hard X-ray observation of the Coma 
cluster region with the IBIS instrument onboard INTEGRAL. In addition to the Coma cluster, the final 20- 
50 keV image contains 12 serendipitous sources with statistical significance > 4<7. We use these data (after 
correcting for expected number of false detections) to extend the extragalactic source counts in the 20-50 keV 
energy band down to a limiting flux of 1.0 x 10~ n erg s _1 cm~ 2 (~ 1 mCrab). This is a more than a factor of 
10 improvement in sensitivity compared to the previous results in this energy band obtained with the HEAO-1 
A4 instrument. The derived source counts are consistent with the Euclidean relation, N(> /) f~ 3 / 2 . A 
large fraction of identified serendipitous sources are low-redshift, z < 0.02 AGNs, mostly of Seyfert 1 type. 
The surface density of hard X-ray sources is (1.4 ± 0.5) x 10~ 2 per square degree above a flux threshold of 
10~" erg s cm~ 2 . These sources directly account for ~ 3% of the cosmic X-ray background in the 20- 
50 keV energy band. Given the low redshift depth of our sample, we expect that similar sources at higher 
redshifts account for a significant fraction of the hard X-ray background. Our field covers only 3% of the sky; 
a systematic analysis of other extragalactic INTEGRAL observations can produce much larger source samples 
and is, therefore, critically important. 

Subject headings: X-rays: general — X-rays: diffuse background — galaxies: active — galaxies: Seyfert 



1. INTRODUCTION 

Most of the energy of the Cosmic X-ray Background (CXB) 
is em itted in the energy band around 30 keV jMarshall et alJ 
1980). However, the exact nature of the source popula- 
tion responsible for the background at these energies is un- 
known. The primary reason is low sensitivity of the previ- 
ous X-ray telescopes operating above 20 keV. Studies of the 
high energy sources are also motivated by the recent work at 
lower energies, 2-10 keV. Most of the X-ray background at 
these energies is reso lved into sources ( Giacconi et al] l2002t 
lAlexander et al.l 2003) but the spectrum of those sources does 
not match that of the CXB at high energies, indicating the 
existence of a population of highly obscured sources, which 
should be detectable more easily above 20 keV (e.g., Worsley 
et al. 2004 and references therein). 

The telescopes onboard INTEGRAL provide a major im- 
provement in sensitivity for X-ray imaging above 20 keV 
(IWinkler et aT1l2003al) . During its first year, INTEGRAL has 
conducted a number of deep pointings to the Galactic cen- 
ter region ( Revnivts ev et alTl2004al) and the Galactic plane 
(Wink ler etaiTl2 003b). as well as to several extragalactic tar- 
gets. One of the deepest extragalactic observations was that 
of the Coma cluster region for a total exposure of 500 ksec. 
The Coma cluster is located very close to the North galactic 
pole and this field is minimally contaminated by the sources 
within our Galaxy which dominate the hard X-ray sky. In ad- 
dition, the image is not "polluted" by a bright target and so it 
is excellent for detection of faint sources. 

In this Paper, we report on the analysis of faint serendip- 
itous sources in the Coma field. Our image reaches a 4(7 
detection threshold of 1 mCrab in the 20-50 keV energy 
band, above which we detected 12 serendipitous sources. 
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Using this sample, we are able to extend the extragalactic 
logN — log S in the hard X-ray band to a flux limit of 1.0 x 
10~ n erg s _1 cm~ 2 , a factor of ~ 10 deeper th an the extra- 
galac tic part of the HEAO- 1 A4 source catalog jLevine et alJ 
119841) . 

2. INTEGRAL OBSERVATION OF THE COMA REGION 

The Coma region was observed by INTEGRAL in 2003 on 
January 29-31 (revolution 36) for 170 ksec and on May 14- 
18 (revolutions 71-72) for 330 ksec. The observations consist 
of 221 shorter pointings which form a 10° x 10° grid around 
the target with a 2° distance between the grid points. The Jan- 
uary and May datasets have different position angles, which 
leads to larger total sky coverage and helps to minimize the 
systematic residuals in the final image. The data quality in 
both series of observations is similar so they can be combined 
in a single ^500 ksec dataset. 

We used the data from IBIS/ISGRI instrument which is 
the most suitable for the imaging surveys in the hard X-ray 
band among the majo r instrument onboard INTEGRAL. IBIS 
dUbertini et alJ l2003) is code d-mask aperture te lescope. Its 
CdTe-based ISGRI detector JLebrun et all 120031) has a high 
sensitivity above 20 keV and has a high spatial resolution. 
The telescope field of view is 28° x 28° (9° x 9° full coded). 
All Coma pointings cover a 40° x 40° region. 

3. DATA REDUCTION AND IMAGE RECONSTRUCTION 

The IBIS data analysis involves special techniques for 
coded-mask image reconstruction, for which we used the soft- 
ware suite developed by one of us (EC). The essential steps 
are outlined below. 

The event energies were calculated following the Off-line 
Scientific Analysis software v3.0 (Goldwurm et al. 2003) us- 
ing the gain table v9 and the event rise time correction v7. Us- 
ing the detector images accumulated in a broad energy band, 
we searched for hot and dead pixels and then screened the data 
to remove these artifacts. This resulted in rejection of several 



2 



KRIVONOS ET AL. 




FIG. 1 . — X-ray image of Coma region in energy band 20-50 keV obtained by IBIS instrument onboard INTEGRAL during two different sets of observation. 
The color-coded intensity is in units of statistical significance for a point source. The markers indicate the locations of sources above a signal-to-noise level of 4. 



percent of the detector area. We then computed raw detector 
images for each pointing location and reconstructed the sky 
images individually. 

The reconstruction starts with rebining the raw detector im- 
ages into a grid with the pixel size equal to 1/3 of the mask 
pixel size. This is close but not exactly equal to the detec- 
tor pixel size. Therefore, the rebining causes a moderate loss 
of spatial resolution but leads to straightforward application 
of the standard coded mask reconstruction algorithms (e.g., 
Fenimore & Cannon 1981, Skinner et al. 1987a). Basically, 
the flux for each sky location is calculated as the total flux in 
the detector pixels which "see" the location through the mask, 
minus the flux in the detector pixels blocked by the mask: 

/= £d-zml£d, (i) 

M=\ M=0 

where / is the source flux, D is the detector image, M = 1 
or corresponds to transparent or opaque mask pixels, re- 
spectively, and BAL w 1 is so called balance matrix. The bal- 
ance matrix accounts for non-uniformity of the detector back- 
ground. For a given location it is calculated as: 

BAL= ZM= 1 D tj (2) 

^M=oDb 

where D\, is the detector image accumulated over large num- 
ber of observations without strong sources in the field of view 



(i.e., the background). Thus t he ex pectation value of / is zero 
for a source-free field (see ea. 11121 . 

The images reconstruction is based on the DLD deconvo- 
lution procedure (see notations in Fenimore & Cannon 1981) 
when the mask pixel corresponds to n x n detector pixels. The 
original detector is treated as n X n independent detectors and 
n x n independent sky images are reconstructed and then com- 
bined into a single image. The point source in such image is 
represented by a n x n square. In our case, this leads to the 
effective Point Spread Function (PSF) being approximately a 
square of 3 x 3 detector pixels or 12' x 12'. 

Periodic structures in the mask produce sky images with 
a number of prominent peaks accompanying every real 
source. An iterative removal procedure was used to eliminate 
"ghosts" associated with the brightest sources (NGC 4151, 
NGC 4388 and Coma). No iterative removal of sources was 
performed for weak sources. The reconstructed images for 
each pointing grid location were co-added in the sky coordi- 
nates. Finally, we computed the statistical uncertainty of flux 
in the reconstructed image. This is a straightforward proce- 
dure because the noise is Poisson in the original data and it 
can be easily propagated through the image reconstruction al- 
gorithm (i.e., data rebining and DLD deconvolution). 

We used the 20-50 keV energy band for detection of the 
serendipitous sources in the Coma field. The choice was mo- 
tivated by the following considerations. The lower boundary 
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of the energy range where IBIS is usefully sensitive is near 
20 keV. The choice of the upper boundary is driven by the 
desire to extend the energy band to the highest energy possi- 
ble without sacrificing the sensitivity. An upper boundary of 
50 keV is a reasonable choice beca use at higher e nergies, the 
IBIS sensitivity rapidly decreases (Lebrunetal. 2003§. We 
checked that at higher energies, there are no sources which 
are not detected in the 20-50 keV band. The reconstructed 
image is shown in Fig. [2 

4. DETECTION OF SOURCES IN THE RECONSTRUCTED IMAGE 

Most of the extragalactic sources are very faint in the hard 
X-ray band. Their detection with INTEGRAL is challenging 
and requires sensitive techniques. Our approach is based on 
the so-called matched filtering technique. A similar method 
applied to the source detection in the ROSAT PSPC images is 
extensively described in Vikhlinin et al. (1995). We refer the 
reader to this paper and present only a brief overview below. 

4.1. Matched Filter 

The matched filter source detection is based on the convolu- 
tion of the original image with the PSF. This method provides 
the most sensitive detection algorithm for finding faint, iso- 
lated sources (e.g., Pratt 1978). The only difference between 
the traditional matched filtering and our approach is to slightly 
modify the filter so that it can provide automatic background 
subtraction. 

The image reconstruction algorithm for coded-aperture 
imaging implies, in principle, that the background in the re- 
constructed image is already subtracted (see eq. [I). How- 
ever, some detector imperfections not fully accounted for by 
the balance matrix and other practical problems might lead 
to residual small-amplitude, large-scale intensity variations. 
The automatic background subtraction is therefore a desirable 
property of the detection algorithm. 

The solution is to build the detection filter, f(r), from two 
components. The positive component in the center matches 
the PSF and provides sensitive detection. The second one 
forms a negative "ring" at larger radii so that the average value 
of f(r) is zero. We decided to implement f(r) at the differ- 
ence between the PSF and a Gaussian of larger width, 

/(r)=P(r)-Aexp(-^, (3) 

where P(r) is the IBIS PSF, and the amplitude satisfies 

A = \ f P{r)rdr. (4) 

a L Jo 

The PSF calibration was obtained from the analysis of Crab 
observations. Crab is essentially a point source for INTE- 
GRAL (the source size smaller than 4' or 1 detector pixel), 
therefore its image can be used as a representation of the 
PSF (we actually used the average of 4 possible reflections 
of the Crab image around the X- and F-axis to reduce statisti- 
cal uncertainties). The derived PSF is essentially a square of 
3x3 pixels (or 12' x 12') surrounded by fainter larger-scale 
wings; its one-dimensional slice is shown in Fig. [2] 

Give n th e PSF, the only free parameter in the filter defined 
by eq. J3I4I > is the Gaussian width, a. The choice of a is driven 
by two conflicting requirements: a) the larger values of a lead 
to better signal-to-noise (SNR) ratios in the convolved images, 
and b) smaller values of a lead to more accurate subtraction of 
small-scale background variations. We chose a = 40' which 
results in the signal-to-noise ratio of 99% of its maximum 
value (that for a — > °°). 
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FIG. 2. — One-dimensional slice through the PSF of our reconstructed 
images after iterative source removal procedure (dotted histogram) and the 
profile of the detection filter (solid histogram). The pixel size is 4'. The PSF 
is the result of coadding of many real images (Crab Nebula) with different 
orientations. The resulting effective PSF is then more round and smooth than 
the ideal PSF (single pointing, no rebinning, perfect detector) which is a flat 
3 by 3 pixels square in the DLD sky. 

4.2. Source Detection, Fluxes, Locations 

The sources are identified as local maxima in the the filter- 
convolved image. To sort out statistical fluctuations, we com- 
puted the convolved noise as 

N(x,y) = ((noise map) 2 ® (/ 2 )) 1/2 . (5) 

This equation is strictly correct only if the individual image 
pixels are statistically independent, which is not fulfilled in 
our case because of rebinning of the detector images (§ 0. 
However, the effect is small because the correlation length of 
the noise in the reconstructed images is smaller than the size 
of the filter. It can be corrected by uniform rescaling of the 
noise map given by eq. l|5} so that the rms variation of the 
ratio of the convolved image and N(x,y) is 1. The required 
correction factors are « 0.9. Note that such adjustment ef- 
fectively corresponds to noise determination directly from the 
data, which is the most trustworthy approach. 

Accepted sources are the maxima in the convolved image 
with the signal-to-noise ratio > 4. This threshold is lower 
than what is often used in wide-area surveys. Our choice was 
motivated by the following considerations. Detection thresh- 
old is the result of compromise between our desire to exclude 
as many false detections as possible (this drives the threshold 
high) while retaining the maximum of real sources in the sam- 
ple (this drives it low). The 4(7 detection threshold results in 
manageable contamination of the sample by false detections, 
4.7 on average or ~ 33% of the total sample, which stil l can 
be reliably subtracted in the final logN — log S analysis (? I6.2I >. 

Source fluxes can be obtained directly from filter-convolved 
images. The peak value in the convolved image is propor- 
tional to total source flux because convolution with the detec- 
tion filter is a linear image transformation. The conversion 
factor between the peak value and flux was obtained using the 
reconstructed Crab images, for which we assumed the Crab 
spectrum 9.7 x E^ 1 phot cm~ 2 s _1 keV -1 l lToor & Seward 
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11974ft . The conversion factor is a function of position because 
of several effects. First, the PSF degradation at large off-axis 
angles changes the peak value in the filter-convolved images. 
Second, off-axis PSF contains weak side-lobes which leads 
to over-subtraction of the source flux by a wavelet-like filter 
with a negative annulus, such as ours. Third, there is s vi- 
gnetting effect most likely caused by absorption in the mask 
holes. These effects were calibrated using a large number of 
Crab observations with the source at different locations within 
the FOV. These observations were reduced identically to the 
Coma field. The derived peak-to-flux coefficient decreases by 
22% from 0° to 5° off-axis and then stays approximately con- 
stant. We used the conversion factor derived at 5° and ignored 
the trend at smaller off-axis angles because it affects only a 
small fraction of the image area. 

The source locations were measured as flux-weighted mean 
coordinates within 6' of the maxima in the filter-convolved 
images. The accuracy of this method was tested using a large 
number of INTEGRAL observations in which there were 
sources with SNR = 5—6. The distribution of coordinate off- 
sets for these sources is well described by a Gaussian with a 
68% uncertainty radius of 4.2'. 

The noise in the reconstructed image increases near the 
edge of the field of view. This region, if included in the survey, 
can pollute the source catalog with high-flux false sources. 
Therefore, we excluded the regions where the rms level of the 
image noise is a factor of 10 greater than that in the center of 
the field of view. This reduced the image area by 28%. 

5. DETECTED SOURCES 

Thirteen hard X-ray sources were detected above the 4a 
threshold in our Coma image, including the target. The 
source locations and observed fluxes in the 20-50 keV energy 
band are given in Tabled We searched for obvious optical 
identifications of our sources using the NASA Extragalactic 
Database as well as by the visual inspection of the images 
from Digitized Sky Survey and the ROSAT All-Sky Survey. 
Seven sources were unambiguously identified with the extra- 
galactic objects (TableQ, most of them classified as Seyfertl 
galaxies at low redshifts, z < 0.02. 

Our brightest sources were observed in hard X-r ays with 
previous observatories (Mac omb & Gehre TsI 1 19991 ). NGC 
4151 wa s observed with HEAO-1 A4 jBaitv et alJJ 19841) . 
SIGMA llMandrou et alll994tlFi"noguenov et alii 1995ft. Bep - 
poSAX JPiroetalJ 119981) . BATSE jParsons et alJ 119981) . 
Seyfert-1 galaxy NGC 5548 was det ected by HEAO-1 A4 
(Rothschil d etalJ 1 1983ft. BeppoSAX ( iNicastro et alJ 12000ft . 
BATSE llLing et alJ 120001). NGC 4388 wa s detected by 
SIGM A ( iLebrun et alJ 11992ft and BATSE ( iBassani et all 
1996). 

Most of positively identified INTEGRAL sources are also 
detected in the soft X-ray band in the ROSAT All-Sky Sur- 
vey. The only exception is the Sy2 galaxy NGC 4388. This 
second-brightest INTEGRAL source is undetectable in the 
ROSAT energy ban d, most likely bec ause of a high intrin- 
sic absorption (Beckman n et alJl20 04). One source (#9) is 
identified with the X-ray bright cataclysmic variable in our 
Galaxy. No identified sources have the same redshift as the 
Coma cluster, the observation target. Therefore, they are not 
located within the cluster or associated large-scale structures 
and we can safely use the entire sample to derive the serendip- 
itous source counts. 

We were unable to unambiguously identify 6 detected 
sources, partly because of the relatively poor positional ac- 



/ [17-50 keV], erg s" 1 cm- 2 
FIG. 3. — Survey area as a function of flux for sources with SNR > 4. 

curacy achievable in the hard X-ray band. We note that some 
fraction of the unidentified sources are likely to be false de- 
tections as discussed below. 

6. LOG N - LOG S FOR DETECTED SOURCES 

6.1. Survey Area 

Since the image noise systematically increases at large off- 
axis distances, our detection threshold of SNR = 4 corre- 
sponds to different fluxes at different locations. Therefore, 
the survey area is a function of flux. This function is com- 
puted easily for the given noise map because the peak values 
in the co nvolv ed image is simply proportional to the source 
flux (see S3. 

Equation (0 was used to obtain the noise map for the filter- 
convolved image. Multiplying this map by 4 and then by the 
flux conversion coefficient gives us the map of the minimal 
detectable flux, / m i n , as a function of position. By count- 
ing the area where / > / m ; n , we obtain the survey area as a 
function of flux, A(f). The results of this computation are 
shown in Fig. [3] The minimum detectable flux at the very 
center of the field of view is 0.9 x 10~ n erg s cm~ 2 (or 
~ 0.9 mCrab). The survey area reaches its geometric limit 
of 1243 square degrees for / > 8.3 x 10~ n erg s cm~ 2 , 
and 50% of this area has the sensitivity better than / = 
1.9 x 10~ 11 erg s _1 cm~ 2 . 

Given the survey area as a function of flux, A(f), the cu- 
mulative source counts (also referred to as the log N — log S 
distribution) can be computed easily as 



N(>f)= ZMftT 



(6) 



fi>f 



However, our source catalog must contain a small number of 
false detections because of the relatively low detection thresh- 
old (SNR = 4). We need to subtract their contribution from 
the logN — log S measurement. 

6.2. Expected Number of False Detections 

The reconstruction algorithm for the coded-aperture im- 
ages leads to the symmetric, Gaussian distribution noise in 
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TABLE 1 

Hard X-ray sources detected in the Coma field 



X-ray ID 


RA 


Dec 


Flux 


SNR 


Optical ID 


Dist. 


Type 


z 




ROSAT 




(J2000) 


(J2000) 


erg s _1 cm -2 






arcmin 






erg s _1 




1 


12 10 33 


+39 24 21 


(3.2±0.1) x 10~ 1() 


44.5 


NGC4151 


0.1 


Syl.5 


0.0033 


7.3 x 10 42 


+ 


2 


12 25 43 


+12 40 42 


(1.1 ±0.1) x 10~ 10 


9.5 


NGC 4388 


1.3 


Sy2 


0.0084 


1.6 x 10 43 




3 


14 17 51 


+25 08 43 


(5.3 ± 1.2J X 10 


5. 1 


JNCjL 5548 


2.0 


Syl 


0.0171 


3.3 X 10 


+ 


4 


12 31 19 


+12 14 01 


(5.2± 1.3) x lO" 11 


4.7 














5 


11 56 47 


+37 00 20 


(3.3±1.0) x lO" 11 


4.0 














6 


13 37 09 


+15 20 00 


(2.9±0.8) x lO" 11 


4.1 














7 


12 18 23 


+29 49 06 


(1.7±0.3) x 10-" 


5.8 


NGC 4253 


0.8 


Syl 


0.0129 


5.9 x 10 42 


+ 


8 


12 59 35 


+27 57 25 


(1.6±0.2) x lO" 11 


7.7 


Coma 1 


3.3 


GClstr 


0.0244 


2.4 x 10 43 


+ 


9 


13 35 12 


+37 12 17 


(1.5±0.4) x 10"" 


4.1 


V* BH CVn 


3.0 


CV 






+ 


10 


13 13 23 


+36 34 20 


(1.6±0.3) x 10"" 


4.6 


NGC 5033 


1.6 


Syl 


0.0029 


2.4 x 10 41 


+ 


11 


12 25 40 


+33 31 13 


(1.4±0.3) x lO" 11 


4.2 


NGC 4395 


2.5 


Syl 


0.0011 


3.0 x 10 40 


+ 


12 


13 10 18 


+20 01 51 


(1.1 ±0.3) x lO" 11 


4.2 














13 


12 48 49 


+33 15 14 


(1.0±0.3) x lO" 11 


4.3 















NOTE. — The uncertainty in the source locations is approximately 4' (90% confidence). The fluxes and luminosities are in the 20-50 keV 
energy band. 



The counts-to-flux conversion for Coma is inaccurate because the source is extended and has the thermal spectrum. INTEGRAL data on 
Coma will be discussed elsewhere. 




Significance 

FIG. 4. — Distribution of the formal statistical significance of the local 
minima in the filter-convolved INTEGRAL image. The solid line shows the 
fit derived from Monte-Carlo simulations of the images with Gaussian noise, 
slightly rescaled to match the observed distribution. 

the resulting image. Therefore, the number of false detections 
(those arising because of noise) above a signal-to-noise ratio 
a can be estimated by counting the local minima with the am- 
plitude below a. 

The small number of the local minima with sufficiently high 
formal significance (SNR < —4) is the main practical diffi- 
culty. We, therefore, adapted the following approach. We 
simulated a large number of images with the Gaussian noise, 
convolved them with the detection filter, and derived the dis- 
tribution of the local maxima as a function of their formal sig- 
nificance, a. The distribution was fit to a third-order polyno- 
mial in the (a,log(dN/do)) coordinates. The analytic model 
derived from simulations was slightly rescaled to match the 
distribution of local minima in the real data (Fig.|4}. 

The rescaled model allows us to predict the number of false 



detections as a function of the statistical significance thresh- 
old. We can also use it to predict the number of false detec- 
tions as a function of limiting flux because we know the im- 
age noise as a function of position. Above our limiting flux, 
1.0 x 10~ n erg s _1 cm~ 2 , we expect 4.7 false detections, and 
1.7 false sources above 2.8 x 10~ n erg s _1 cm~ 2 . The ex- 
pected average number of false detections is less than 0. 1 at 
fluxes / > 10~ 10 erg s _1 cm~ 2 . We can also compute the cor- 
responding logN — log S function for false sources using the 
survey sky coverage, identically to what is done with the real 
sources, 

AW(>/) = Lf ^da, (7) 

7-/0..,., A: Oil;., I 

where n, ; is the image noise in the pixel Cmm = 

max(4,//n I j), A(f) is the sky coverage as a function of flux, 
dN I da is the number of false detections in one pixel as a 
function of the formal statistical significance, and the sum is 
over all image pixels. 

This procedure for estimating the number of false detec- 
tions has been checked using shorter, 64 ksec, pieces of 
the Coma observation. Those sources detected in the short- 
exposure image and absent in the total 500 ksec image, are 
false. We verified that the logAf — log S for false sources pre- 
dicted by eq. for short exposures is in excellent agreement 
with that observed. 

6.3. Results 

Figure [5] shows the logN — log S distribution for the 
serendipitous sources in the Coma field, corrected for the 
survey sky coverage as a function of flux, and subtracted 
contribution of false detections. Our results can be com- 
pared directly only with the previous measurements with the 
HEAO-1 A4 instrument which operated in the overlapping 
energy band, 25-40 keV JLevine et alll984l) . 

The logN — log S for the HEAO-1 A4 source catalog in the 
extragalactic sky (|fe| > 30°) is shown by the solid histogram 

in Fig. [5] A single Euclidean function, N(> /) « /~ 3,/2 , can 
fit both datasets. This is not surprising given that our sources 
are at low redshifts. For the fixed slope of the logiV — log S 
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FIG. 5. — The logiV — logS distribution for the serendipitous hard X-ray 
sources in the Coma field. The histogram shows the source counts derived 
from the HEAO-1 A4 extragalactic sample (\b\ > 30°). The HEAO-1 A4 
fluxes were converted to our energy band from the nominal 25^10 keV band 
assuming a power-law spectrum with T =2. The dashed lines shows aAI« 
/~ 3 / 2 fit to our dataset. 

distribution at —3/2, we derive from our data only the sur- 
face density of the extragalactic hard X-ray sources N = 
(1.4 ± 0.5) x 10~ 2 per square degree above a limiting flux of 
10 -11 erg s _1 cm~ 2 in the 20-50 keV energy band. 

All identified INTEGRAL sources have low redshifts, z < 
0.02. Therefore, the depth of our survey is small and it can 
be affected by the nearby large-scale structures. It would be 
extremely important to extend such measurements to a larger 
number of fields. However, it is interesting to note that, thanks 
to a greater sensitivity of INTEGRAL, the volume covered 
by our observation is larger than the volume covered by the 
extragalactic (\b\ > 30deg) portion of the HEAO-1 A4 all-sky 
survey. 

The all- sky catalog from the RXTE slew survey 
jRevnivtsev et alJ l2004bl) reaches a similar depth as our 
survey, although at a lower energy band of 8-20 keV. We 
can compare the number densities of sources by converting 
the RXTE fluxes to our energy band assuming a power law 



spectrum with T = 1.6, the typical photon index for the RXTE 
AGNs (Fig. 8 in Revnivtsev et al.). The RXTE \ogN - log S 
predicts the source density of (0.9 ±0.1) x 10~ 2 per square 
degree above the flux limit of our survey; this is somewhat 
lower than, but within the uncertainties of, our measurement. 



7. DISCUSSION AND CONCLUSIONS 

We have presented the analysis of the deepest hard X-ray 
image of the extragalactic sky obtained to date. We detected 
12 serendipitous sources in the 40° x 40° field centered at the 
Coma cluster. Most of the detected sources appear to be the 
Syl galaxies at low redshifts, z < 0.02. 

The logN — log S distribution derived from a combination 
of our survey and extragalactic sample of the HEAO-1 A4 
source catalog follows the Euclidean function, N(> /) °< 
y-3/2 rp ng observed normalization of the logN — log S, 
(1.4 ±0.5) x 10~ 2 per square degree above a limiting flux 
of 10~ n erg s _1 cm~ 2 , corresponds to the integrated flux of 
4.3 x 10~ 13 erg s _1 cm~ 2 deg~ 2 , or 3% of the total intensity 
of the hard X-ray bac kground in the 20-50 keV energy band 
(Marshall eTallFjSOl) . 

The redshift depth of our source catalog appears to be be- 
low z = 0.02. Therefore, it is reasonable to expect that the 
Euclidean source counts will extend to much fainter fluxes 
than our sensitivity limit. For example, the depth of z < 0.2 
where the space curvature of AGN evolution effects are still 
small, will correspond to approximately a factor of 100 lower 
flux. Extrapolation of our logiV — logS by a factor of 100 to- 
wards fainter fluxes will account for 30% of the total X-ray 
background in the 20-50 keV energy band. In short, we start 
to uncover the source population responsible for a significant 
fraction, if not most, of the hard X-ray background. 
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